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ABSTRACT 



A certain representation for the Heisenberg algebra in finite- 
difference operators is established. The Lie-algebraic procedure of 
discretization of differential equations with isospectral property 
is proposed. Using s/ 2 -algebra based approach, (quasi)-exactly- 
solvable finite-difference equations are described. It is shown that 
the operators having the Hahn, Charlier and Meixner polynomi- 
als as the eigenfunctions are reproduced in present approach as 
some particular cases. A discrete version of the classical orthog- 
onal polynomials (like Hermite, Laguerre, Legendre and Jacobi 
ones) is introduced. 
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Discretization is one of the most powerful tools of solving continuous the- 
ories. It appears in different forms in various physics sciences leading to 
discrete versions of differential equations of classical mechanics, lattice field 
theories etc. However, once the discretization is chosen as an approach to a 
problem we meet a hard problem of ambiguity - there exist infinitely-many 
different ways of discretization having the same continuous limit. The goal 
of present paper is to demonstrate a certain discretization scheme for differ- 
ential equations with exceptional properties - eigenvalues remain unchanged 
(isospectrality), eigenfunctions are modified in a simple manner etc. 

1. Two operators a and b obeying the commutation relation 

[a, b] = ab — ba = 1, (1) 

with the identity operator in the r.h.s. define the Heisenberg algebra. The 
Heisenberg algebra plays the central role in many branches of theoretical and 
mathematical physics. A standard representation of (1) in the action on the 
real line is the coordinate-momentum one: 

= (2) 

Our goal is to build the one-parametric representation of (1) on the line in 
terms of finite-difference operators having (2) as a limiting case. 

Let us introduce the finite-difference operators T>± possessing the property 
of the translation covariance 

v+s(x) = /(* + *) -/(,) s (j^i) 



and 



f(x)-f(x-S) _ (1-e 



-<s-f ■ 

ax 



V_f{x) = J " ' y = ± '-f{x) (4) 

and V + — > D_, once 5 — * —5. It is worth noting that 

V+-V- = 5V_V + (5) 

Now define the so-called quasi-monomial 

= x(x-5)(x-2S) ...(x-n8) = , (6) 
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which in the limit 5—^0 becomes the ordinary monomial x n+1 . The operators 
T2± (3)-(4) act on quasi-monomials as follows 

V+x^ = nx {n - 1] , x{l - 5V^)x (n) = x in+1 ^ . (7) 



So the operator V + plays the same role in action on the quasi-monomial 
x^ as the operator of differentiation in action on monomial x n : it converts 
nth quasi- monomial to (n— l)th with the coefficient of proportionality equals 
to n. The operator x(l — 8T>_) acts similarly to the operator of multiplication 
x on monomials: it converts nth quasi-monomial to (n+l)th quasi- monomial. 
In the limit 5^0 those operators become and x , respectively. Straight- 
forward calculation shows that the commutator of the operators T> + and 
x(l — 5T>_) equals to one ! So the operators 

a = V + , 

b = x(l-SV_) (8) 

form one-parametric representation of the Heisenberg algebra (I). It is easy 
to show that (8) belongs to the point canonical transformations. 

One of possible interpretations of (8) is to identify the finite-difference op- 
erator a with momentum in discretized form, while b becomes the conjugated 
coordinate (in a complete agreement to the limit 5 — > 0). 

2. Consider the eigenvalue problem for some linear differential operator 

L ^^ x \ = Xl P( x ) ( 9 ) 

having polynomial eigenf unctions. Now let us replace in the equation (9) 
the derivative ^ and x by the operators a and b, respectively, obeying the 
commutation relation (f). Then we arrive at 

L[a,b](p(b) = \<p(b) (10) 

In order to make sense to (10), one should add the definition of the vacuum 
|0 >: 

o|0 > = . (11) 
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Then the equation (10) has a meaning of the operator eigenvalue problem 
in the Fock space with the vacuum (11). It is evident that the eigenvalue 
problem (10)— (11) has the same eigenvalues and the same eigenf unctions 
(with replacement of x —> b) as the original problem (9). 

So taking different realizations of the algebra (1) and considering a cer- 
tain element of the universal enveloping algebra of the algebra (1), we get 
different eigenvalue problems with the property of isospectrality. Below we 
will consider isospectral problems emerging from the representations (2) and 
(8). We name such a way of discretization the Lie-algebraic discretization. 

Now we have to specify the operators L in (10) having polynomial eigen- 
functions. In [|IJ it was proven that L has a certain amount of polynomial 
eigenfunctions if and only if, L is the superposition of the element of the uni- 
versal enveloping algebra of the s/2-algebra taken in the finite-dimensional 
irreducible representation 

= b 2 a — nb 
n 

J° n = ba-- (12) 
J" = a 

where n is a non-negative integer, and the operator B(b)a n+1 , where B(b) is 
any operator of b. The dimension of this representation is equal to (n+1) and 
hence (n + 1) eigenfunctions have a form of polynomial of degree n. These 
operators L are named quasi-exactly-solvable. Moreover, one can prove that 
L possesses infinite-many polynomial eigenfunctions if and only if L is a 
polynomial in generators J° = J®, J~ = Jq only Those operators L are 
named exactly-solvable. 

3. In order to exploit the representation (8) let us firstly define the vacuum 
|0 >. The condition (11) in the explicit form is 

f(x + 5) = f(x) . (13) 

Any periodic function with the period 5 is the solution of this equation, 
however, without the loss of generality we can make the choice 

f{x) = 1 . (14) 
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With this vacuum, it is easy to see that 

b n \0 >= [x(l - 5£L)] n |0 >= x {k) (15) 

where the relation 

[ xe - 5 £} n = x^e~ n5 £ . 

was used for derivation (15). So, starting from (10)— (11), we obtain the 
eigenvalue problem for the finite-difference operator 

L[V+,x(l -8VJ)](p{x) = \<p(x) (16) 

where the solutions (p{x) are related to the solutions 

(p(x) of (9). To clarify this relationship, let us assume that 

ip{x) = ]T a kX k (17) 

is a certain solution of the equation (9). The transition from (9) to (10) 
implies the replacement of x by b. Then taking into account (15) we come to 
the conclusion that each monomial in (17) should be replaced by the quasi- 
monomial. Finally, we get that the corresponding solution of the equation 
(16) has the form 

^x) = J2a k x^ (18) 

The equation (16) has the same eigenvalues as the original equation (9), 
while the eigenfunctions are modified by replacing the monomials by quasi- 
monomials. It is evident, that such a procedure preserves the degree of 
polynomial. 

Now let us proceed to concrete realization of idea of the Lie-algebraic 
discretization for differential equations with polynomial solutions. It is quite 
obvious that there exist two different approaches: 

• In the first approach, we start from a certain differential equation 
with known spectra and apply above-described Lie-algebraic discretiza- 
tion. Generically, a second-order differential equation becomes a finite- 
difference equation relating the unknown function in more than three 
points or, in other words, more than the three-point finite-difference 
equations. 
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• The second approach is different: we look for the differential equa- 
tions leading to the three-point finite-difference equations under the 
Lie-algebraic discretization. In other words, what are the elements of 
the universal enveloping s/ 2 -algebra in the representation (8), (12) pro- 
ducing the three-point finite-difference equations. 

4. The first approach. In |]J (see also 0) it is shown that the most general, 
second-order exactly-solvable differential operator on the real line possessing 
infinitely-many polynomial eigenfunctions has the form of the hypergeomet- 
rical operator 

E 2 (±x) = -Q 2 ( X )*L + Ql ( x )*.+Q (19) 

where Qk{x) are polynomials of the kth order with arbitrary coefficients: 

Q 2 = a x 2 + a\x + a 2 , Q\ = b x + bi, Q = c . 

The operator (19) is equivalent to the operator appearing from the most gen- 
eral quadratic polynomial in the generators J°, J~ (see (12)) defined through 
the representation (2) [|IJ . Performing the above-mentioned procedure of the 
Lie algebraization: (9) — > (10) — > (16), we arrive at the following exactly- 
solvable finite-difference operator isospectral with (19) 

E 2 = E 2 [V + ,x(l-5V_)] = 

-a 2 e 25 ^ + [-a x x + (2a 2 + bi)}e 5 ^ + 
[— a x(x - 6) + (2ai + b )x - (a 2 + b\)] + 
[2a x(x -5) - (at + b )x}e' s ^ - a x(x - 5)e" 2<5 ^, (20) 
where on = ^i,bi = h. The corresponding eigenvalue problem has the form 

-a 2 ip(x + 25) + {-a x x + (2a 2 + l$\<p{x + #)+ 

[— a x(x - 5) + (2ai + b )x - («2 + bi)](p(x)+ 
[2a x(x — 5) — (5i + b )x](p(x — 5) — a x(x — 5)<p(x — 26) = X(f(x), (21) 
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It is worth noting that in contrary to a naive expectation, although it is 
started from the second-order differential operator, the finite-difference op- 
erator connects to the function in five different points: (x + 25, x + 5, x, x — 
S,x-2S). 

Now taking, for instance, a concrete operator (19) with the Hermite poly- 
nomials H k (x) as the eigenfunctions (a = a\ = b\ = 0, a 2 = — 1, b = —2): 

d d 2 ^ d 

dx : dx 2 dx' 

or, equivalently, 

h = J~J~ - 2J° 

and performing the Lie algebraization (taking into account (5)), we arrive at 
the isospectral finite-difference operator 

h = V\-2xV_ (23) 

The corresponding eigenvalue problem has the form 

1 2 1 1 2x 

—<^(x + 25) - —<p(x + 5)- ~{2x - -)y(x) + -j^ x _ 5 ) = A ^0)' 

with the eigenvalues = 2k and the eigenfunctions H^(x) = Y?iZoCi x \ 
where k = 0, 1, 2 . . . and q are the standard Hermite polynomials coefficients. 
We name these polynomials H^(x) the discrete Hermite polynomials. In 
similar way one can construct the discrete Laguerre, Legendre, Jacobi poly- 
nomials. The results will be presented elsewhere. 

Analogously, one can take the general second-order quasi-exactly-solvable 
operator 

T 2 {^-,x) = -P,(x)^ + P 3 (x)^ + P 2 (x) (24) 

where Pk{x) are polynomials of the kth order depending on ten free param- 
eters one of which is a non-negative integer, n (see (12); for more details 
see 0). It leads to the following isospectral, quasi-exactly-solvable, finite- 
difference operator 
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f 2 = T 2 [V+,x(l-5V-.)} = 

a 2 e 2S ^ + P h i(x)e 5 ^ + 
P 0;2 (x) + xP_ 1)2 (x)e- 5 s + x^P^ 2t2 {x)e- 25 ^ + 

xWp„ 3A (x)e- 3S £ + a_ 4 x (4) e- 45 ^, (25) 

where a's are parameters and Pj^{x) are polynomials of the fcth order. It is 
worth emphasizing that the quasi-exactly-solvable finite-difference equation 
corresponding to (25) relates unknown function to seven different points : 
(x + 25, x + 5, x, x — 5, x — 25, x — 35, x — 45). 

5. The second approach. Standard second-order finite-difference equation 
relates an unknown function at three points and has the form 

A(x)(p(x + 5)- B(x)(p(x) + C(x)(p(x -5) = \ip(x), (26) 

where A(x), B(x), C(x) are some functions. One can pose a natural problem: 
what are the most general coefficient functions A(x), B(x), C(x) for which the 
equation (26) admits infinitely-many polynomial eigenf unctions ? Basically, 
the answer is presented in ]lj: any operator with above property can be 
represented as a polynomial in the generators : 

J°=§(l- e -^) 
o 

J- = I(e 5 ^_i) (27) 
5 

One can show that the most general polynomial in the generators (27) leading 
to (26) is 

E = A 1 J° J°( J- + -)+ A 2 J°,r + A 3 J° + A 4 J~ + A 5 (28) 
5 

and in explicit form 
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,Aa Ao A\ q_ x d 

[ T + ^ x + ^ x]e + 

J-\A Ai J_ ( Al ^0 A * A3 \ Al 21 

+[A5 -T + ( ^ + V~T )x -V x] 

r , A-\ Ao A"t . A-\ o, c d , . 

+ ["(^-^ + f )X+ i x]e (29) 
The spectral problem corresponding to the operator (29) is given by 

+ - ^ + t' 1 + ^/(z - *) = A /w ■ < 3 °) 

In general, this spectral problem has the Hahn polynomials (x, N) as 
the eigenf unctions (therein we follow the notations of ||). Namely, these 
polynomials appear, if 5 = — 1, A$ = and 

A = -1, A 2 = N — (3 — 2, A 3 = -a - (3 - 1, A 4 = {(3 + 1)(N - 1) . 

Besides that, if 

A 1 = 1, A 2 = 2 - 2N - u, A 3 = 1 - 2N - ^ - p, A 4 = (N + u - 1)(N - 1) 

the so-called analytically-continued Hanh polynomials h^' v \x,N) appear, 
where k = 0, 1, 2 . . .. 

So the equation (30) corresponds to the most general exactly-solvable 
finite-difference problem, while the operator (28) is the most general element 
of the universal enveloping s^-algebra leading to (26). Hence the Hahn 
polynomials are related to the finite-dimensional representations of a certain 
cubic element of the universal enveloping s^-algebra (for a general discussion 
see§). 

Taking in (30) S = 1, A$ = and putting 

A x = 0, A 2 = -/i, A 3 = \i - 1, At = 7/i , 
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we reproduce the equation having the Meixner polynomials as the eigenfunc- 
tions. Furthermore, if 

A 1 = 0,A 2 = 0, A 3 = -1,A 4 = n , 

the equation (30) corresponds to the equation with the Charlier polynomials 
as the eigenf unctions (for the definition of the Meixner and Charlier polyno- 
mials see e.g. ||). For a certain particular choice of the parameters, one can 
reproduce the equations having Tschebyschov and Krawtchouk polynomials 
as the solutions. 

Among the equations (26) there also exists quasi-exactly-solvable equa- 
tions possessing a finite amount of polynomial eigenf unctions. All those 
equations are classified via the element of the universal enveloping s/2-algebra 
taken in the representation (8), (12) 

T = A f (J+ + M£j£)+A 1 J^ (31) 

(cf. (28)), where ^4's are free parameters. 

In closing the authors would like to express their gratitude to the Instituto 
de Fisica, UNAM for kind hospitality extended to them. This work was 
supported in part by the research CONACyT grant. 
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